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1. INTRODUCTION 

During January 1981 the U.S. Air Force (USAF) Flight Test Center 
(AFFTC) at Edwards Air Force Base, California, conducted a test program at 
the NASA Lewis Research Center (LeRC) Icing Research Tunnel ( IRT) aimed at 
documenting and improving the performance of the USAF A-10 inlet Icing 
Protection System. While inlet testing was the primary objective, the 
testing period and schedule provided a unique opportunity for additional 
icing instrumentation studies. These additional studies were focused on 
comparisons between modern state-ofthe-art droplet/icing measurement sys- 
tems and those systems and techniques already in use for some years. In 
addition, documentation and investigation of the characteristics of the 
IRT itself could be accomplished with both types of systems. 

These studies were beyond the scope of interest to the USAF A-10 
group. However, they were of vital interest to NASA, FAA, the U.S. Army 
and Air Force and indeed to the entire icing and cloud physics community. 
In August 1981, MR I was awarded a contract (NAS3-22760) from NASA LeRC 
to perform the analyses of the data obtained during USAF test program in 
January. The primary objectives of the contract were to provide the 
followi ng: 

• Descriptions of the "modern" icing instrumentation 
i nstal led i n the IRT. 

• Data listings for each tunnel test condition. 

• Analyses of the data including information investi- 
gations into: 

1. The repeatability of the tunnel and tunnel 
uni formi ty, 

2. Instrument comparisons with other instruments 
and comparison of tunnel droplet spectra 
with natural cloud conditions, and 

3. Velocity effects on spectrometers and temper- 
ature effects on accretion instruments. 


1 


This final report serves to accommodate these objectives in fulfill- 
ment of the report requirements of the contract. Section 2 presents descrip- 
tions of the instruments in addition to the instrumentation setup in the IRT. 
Documentation of the data reduction procedures for the data listings which 
are presented in a separate data volume is given in Section 3. Data recovery 
information and problems encountered during testing are also discussed. The 
results of the analyses are presented in Section 4 with conclusions and 
recommendations following in Section 5. 


2. IRT INSTRUMENTATION AND DESCRIPTIONS 

Figure 2-1 presents a schematic of the instrumented IRT. In addition 
to the instruments for recording of tunnel air speed, air temperature, and 
dew point temperature, the test section was equipped with a specially fabric 
cated test stand for mounting of the cloud droplet instrumentation. The test 
stand was equipped with a Johnson and Williams hot wire deyice and a Leigh 
ice detection meter and was capable of mounting two of the Particle Measure- 
ment Systems, Inc. (PMS) cloud particle spectrometers as shown in Figure 2-2. 
The PMS probes tested included four Axially Scattering Spectrometer Probes 
(ASSP), two Forward Scattering Spectrometer Probes (F$SP), and two Cloud 
Particle Spectrometers (CPS). Probes belong to the following organizations 
were used in the IRT: 

• U.S. Air Force Geophysical Laboratory (AFGL), 

• Meteorology Research, Inc. (MRI). 

• Particle Measurement Systems, Inc, (PMS) (leased by MRI). 

• Leigh Instruments, Limited. 

The FSSP's and CPS's were Interchanged at the location shown for the CPS in 
Figure 2-2. The ASSP's were always mounted as shown in the figure. Data 
were recorded unto digital tape at a rate of once per second. 

Since the heaters provided with the PMS probes were not able to de-ice 
the probes adequately at the cold temperatures, the PMS probes were wrapped 
with steam tubes. Special tests revealed that the steam tubes did not alter 
the air flow through the probes' sample tubes. Also, for given test condi- 
tions, near identical estimates of the cloud properties were determined from 
the measurements with and without the steam tube wrappings as will be dis- 
cussed in a later section. 

The following subsections present descriptions of the probes installed 
on the test stand. 
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Figure 2-1. TOP VIEW OF IRT SHOWING INSTRUMENT LOCATIONS 
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2.1 Axially Scattering Spectrometer Probe (ASSP-100) and Forward Scattering 

Spect rometer Probe (FSSP~-iOQ) 

The Particle Measurement Systems, Inc. (PMS) ASSP-100 is a droplet siz- 
ing device developed in the early 1970 ' s. Droplet sizes are estimated on a 
particle-by-particle basis from measurements of the amount of forward scattered 
light as droplets pass through a focused laser beam. Figure 2-3 presents a 
block diagram of the optical system. Particles intercepting the laser beam 
scatter light through the optics. The scattered light is received at a beam- 
splitting prism where one half of the beam is deflected 90° and focused onto 
a detector (annular) which is masked at the center while the other half is 
focused onto an unmasked (or signal) detector. 

Since the intensity of light scattered from a given size droplet is a 
function of its position within the focused beam, it is necessary to establish 
if the particle is contained within the desired depth-of -field and if the par- 
ticle passes through the central part of the beam where the intensity is 
uniform. The depth-of-field is defined as that region along the beam axis 
where the signal received by the unmasked detector exceeds that of the masked 
detector by a factor of 2.5. A droplet is rejected if the pulse pair compar- 
isons indicate a smaller value. The transient time (or pulse width) of the 
particle as it passes through the laser beam is analyzed to determine its 
chordal transect through the beam. Particle transit times are measured and 
compared with a stored running mean of approximately 1,000 accepted pulse 
widths. Pulses exhibiting transient times greater than the average are 
accepted. This edge effect discrimination is performed to eliminate under- 
sizing of particles as they pass through the outer and weaker intensity 
regions of the beam and is based upon the decrease in chordal transit dis- 
tances from the center through the edge of the beam. 

For a particle accepted as being within the depth-of-field and in the 
central region of the beam, the peak intensity is measured by a pulse height 
analyzer to determine particle size. Since Mie theory for spherical parti- 
cles indicates that scattering intensity is near proportional to the square 
of particle size, the signal from the pulse height analyzer is gated to pro- 
vide a count in one of 15 equally spaced size channels by a series of voltage 
comparators having square weighted voltage thresholds. 
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Figure 2-3. AXIAL SCATTERING SPECTROMETER PROBE OPTICAL SYSTEM DIAGRAM 





Each particle sized and counted by the ASSP causes a latch to be set 
and inhibits the measurement of other particles for a period of 4 /us. This 
finite time period insures sufficient time for quantification and transfer of 
information without interference from subsequent particles entering the samp- 
ling aperature* 

The FSSP was designed and built as an improved successor to the ASSP. 
The modifications included refinements to its optical and electronic sub- 
systems rather than in basic operating principles. As compared to the ASSP, 
the FSSP analyzes light scattered at wider angles of 4-18° as opposed to the 
7-16° in the ASSP. The FSSP also generally has a shallower depth of field 
which results in sharper size discrimination and a smaller sampling area. 

The physical construction of the FSSP is also different as shown in the block 
diagram of the optical system in Figure 2-4. 

Both probes are capable of sizing particles from 0.5 to 47 microns 
diameter in four size ranges {2=47 jum, 2-32 ftm, 1-16 fim, and 0.5-8 jum) and 
resolves them into 15 equally spaced size intervals. They are also normally 
designed for airborne sampling velocities of 20 to 125 m/sec (~250 kts) and 
can be modified for use at lower or higher velocities. For measurements con- 
ducted at higher speeds, pulse height distortion leads to under estimation 
of particle size. The characteristics of the ASSP and FSSP are summarized 
in Tables H«‘ and 2-2. As an option* the total number of pulses or strobes 
sensed through the beam width and within the depth-of-field can be recorded. 
Hence, an acceptance ratio (or efficiency) based on the total particle counts 
in the 15 size channels and total strobes provides a more accurate value for 
the effective beam diameter. Early versions of the ASSP did not have this 
velocity reject feature and, therefore, assumed an effective diameter of 
approximately one half the beam width. 

Another option on probe activity serves to provide information on 
possible errors in the counting efficiency of the probes. Since particles 
which are sensed as entering the beam (outside as well as within the depth 
of field) require finite time periods in the accsptance/rejection logic, they 
are also counted and processed to provide outputs on probe activity from 
which the "busyness” of the probe is inferred. High (as opposed to low) 
activity indicates that a 'large (small) portion of the particles entering 
within the effective sampling area (effective beam width x depth of field) 
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Figure 2-4. FORWARD SCATTERING SPECTROMETER PROSE OPTICAL SYSTEM DIAGRAM 





Table 2-1 


AXIALLY SCATTERING SPECTROMETER PROBE CHARACTERISTICS 


Number of size channels 
Size range 
Size resolution 
Maximum particle rate 
Coincidence errors 
Maximum particle velocity 
Power 

Envi ronmental 

Mechanical 

Weight 


15 

3 to 45 microns 
3 microns 
100 k Hz 

Less than 10% with Concentrations of 1G 3 cm” 3 
125 m/sec 

115v, 60 Hz, 1.0 amp plus required de-ice 
power (300 W using 28 v heaters) 

Temperature: -40°C to +40°C 

Altitude: 0 to 40,000 ft msl 

Probe: 24-inch long airfoil 

Probe: approximately 6 lb 


Table 2-2 

FORWARD SCATTERING SPECTROMETER PROBE CHARACTERISTICS 


Number of size channels 
Size range 
Size resolution 
Maximum particle rate 
Coincidence errors 
Maximum particle velocity 
Power 

Environmental 

Mechanical 

Weight 


15 

3 to 45 microns 
3 microns 
100 k Hz 

Less than 10% with concentrations of 10 3 cm” 3 
125 m sec”* 

115v, 60 Hz, 1.0 a plus required de-ice power 
(300 W using 28 v heaters) 

Temperature: -40°C to +40°C 

Altitude: 0 to 40,000 ft msl 

Probe: 34-inch long, 7-inch diameter cylinder 

Probe: approximately 45 lb 
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would not be processed since a large (small) fraction of a given time period 
is lost in rejecting particles. Laboratory experiments at PMS show that the 
counting efficiency (probe count/true count) decreases exponent telly from a 
value of 1 starting at an activity of about 40 percent to approximately 0.6 
as activity increases to 80 percent. Hence bulk parameters such as liquid 
water content (LWC) and total number concentrations are usually underestima- 
ted at the higher activities. This type of coincidence error is thought not 
to affect the resulting size distribution shape or form from which medium 
volume diameter (MVD) is estimated since it applies equally to all particles. 
Although PMS provides empirical correction factors from their laboratory 
tests, it is felt that since scattering intensity is dependent upon particle 
size, the corrections would not necessarily apply to clouds having different 
particle size distributions. This contention is partially verified in Sec- 
tion 4 where LWC estimates from probes having this option are compared. 

2.2 Optical Array Cloud Particle Spectrometer Probe (0AP-200X) 

The Cloud Particle Spectrometer counts and sizes particles from 20 to 
300 fx m by imaging the shadow of each particle on a linear array of photodiodes 
as shown in Figure 2-5. A helium-neon laser (6238ft) is used with condensing 
and imaging optics to focus the shadow of a particle lying within the sample 
volume. The beam is vertical between sampling arms, and the focal plane is 
centered between the two sampling arms. The particle may shadow one or more 
of the photodiode elements. Each active element of the photodiode array con- 
tains its own amplifying and logic circuit (PC). If the shadowing of an ele- 
ment causes at least a 50 percent reduction in light level to PC (this is an 
electronic truncation of depth of field and minimum size), a flip-flop circuit 
is triggered. Two end photodiode elements are used in logic to reject parti- 
cle shadows extending over the end of the active part of the photodiode array. 

Information from each PC of the array is acquired by an accumulator. 

The particle is sized in the BMS by accumulating flip-flop information from 
the PC, and converting this inco a particle size. The unit particle size 
results from a determination of the number of elements set, tne size of each 
array element, and the magnification of the optical system. Each second of 
particle information collected by the BMS is dumped in proper form for 
recordi ng. 
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Figure 2-5. CLOUD PARTICLE SPECTROMETER PROBE OPTICAL SYSTEM DIAGRAM 


These probes are equipped to size into 15 equal size channels. There 
are two basic size ranges currently in use. The cloud probe range is 20 to 
300 /im, and the precipitation size range generally over 200 to 4500 /;m. In 
the precipitation size range, some models size 200 to 3000 /im, and others 300 
to 4500 Mm, each divided into 15 equal size channels. Detailed specifications 
of these probes taken from PMS handbooks appear in Table 2-3. 

Table 2-3 

CLOUD PARTICLE SPECTROMETER CHARACTERISTICS 


Number of array elements 24 

Number of size channels 15 


Size range (basis system) 
Maximum particle velocity 
Maximum aperture 
Power 

Dimensions 

Weight 

Envi ronmental 


20 to 300 microns 
250 knots 
6.5 cm 

115 v, 50 to 400 or 60 Hz, 60 W 
De-ice: 28v DC, 70 W 

Cylinder: 18 in (length); 6-1/2 in (diameter); 

Optical Extensions (two): 20 in (length); 

1 in (diameter) 

45 lb 

Temperature: -60°C to +50°C 

Humidity: 10% to 100% 

Altitude: 0 to 40,000 ft msl 


Of importance to users of these probes is their sampling volumes. 
Knollenberg (1970) used glass beads and opaque discs of various diameters 
placed on glass slides to determine their shadow size and light intensity 
distribution in a parallel plane as a function of distance. Using coherent 
illumination, he found that at close distances, the particle was in sharp 
focus and all of the decrease in light intensity was measured to be in a 
distance equal to the particle diameter. One would expect this from basic 
optics (Stone, 1963) and for non-spherical solid articles one would expect 
that at close distances all of the decrease in light intensity would be in 
an area equivalent to the cross sectional area of the particle. At distances 
far from the particle diameter, the decrease in intensity was found to be 
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distributed over a much wider distance than the particle diameter. The 
depth of field DF for the measured shadow size to be within ± 10 percent of 
the actual particle size, which corresponded to a 40 percent reduction in 
light intensity, was found in non-dimensional units to be DF * ± 3. The 
detection threshold of a 50 percent reduction in light level before the cir- 
cuits are triggered was derived from this criteria. The depth of field In 
dimensional units as a function of particle radius R and illumination wave- 
length X is therefore, DF = + 3R 2 /x. For a He-Ne laser where X « 0.6328 pm, 
for a 200 pm diameter particle DF «= 9,4 cm, and for a 1000 Xm diameter parti- 
cle, DF B 948 cm. 

The sampling volume is the product of the depth of field, the effective 
array width (which Is related to the width of the array minus the diameter of 
the particle) and the true air speed. If the depth of field for a given par- 
ticle diameter is wider than the >eparation between sampling arms, then the 
depth of field is equal to the arm width of 6.1 cm. This is true for the 
drops having diameters larger than 160 pm. 

2.3 Leigh Ice Detector Set 

The Leigh Ice Detector Set (IDS) aspirator features an annular ejector 
nozzle which whan powered by compressed air is used to entrain ambient air 
and thereby induce high-velocity air flow over the ice-collecting probe 
(Figure 2-6). The air flow is maintained relatively constant despite varia- 
tions in the freestream velocity vector. Thus, the IDS is designed to con- 
tinuously sample the air for supercooled liquid water. 

The ice-collecting probe which is located crosswise in the ejector 
duct takes the form of a smal 1 -diameter thin-wall tube (Figure 2-7). It has 
an ice-collecting efficiency close to unity as calibrated in the Canadian 
National Research Council's Icing tunnel. The thin-wall tube has sufficient 
resistance to serve as its own heater with no internal filament required. 

The thermal inertia of the probe is therefore very low. In response to 
heat generated by probe de-icing, the ice interface at the probe is melted 
and the ice is shed out the tailpipe of the duct. The photodetector 
senses restoration of the beam to full strength at shedding and turns 
the probe de-ice power off. Continued cyclic icing and de-icing occurs 
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until no further ice is encountered. It is notable that the presence of 
supercooled droplets is determined from icing of the supercooled water and 
not just the atmospheric conditions under which ice can but might not form. 

Crosswise in the duct and at a 45° angle to the probe axis is located 
a solid-state soft infra-red emitter and a photo transistor (Figure 2-7). 

The oblique light path ensures that a large ice catchment area on the sensor 
probe face is monitored for ice growth thereby accommodating minor variations 
in probe ice formations. Once a given thickness of ice is accreted on the 
probe, the IR beam is occluded and logic circuits within the unit turn on 
electric current to de-ice thr probe. 

The sawtooth output signal from the photo transistor resulting from 
the cyclic buildup of ice on the probe is then electronically processed to 
provide an initial ice warning signal and a computation of rate of formation 
of ice in terms of ambient liquid water content (LWC). 

The ice detector unit provides a voltage signal that is proportional 
to the liquid water content being measured. This signal is used to drive a 
standard aircraft type instrument scaled in liquid water content (LWC). The 
signal can also be recorded for data gathering purposes. 

2.4 Johnson-Wi 1 liams (J-W) Liquid Water Content Indicator 

The Johnson-Williams (J-W) liquid water content indicator is a hot- 
wire aircraft instrument designed to measure the liquid water content (LWC) 
of droplets having diameters of less than about 40 /im and is an adaptation 
of the hot-wire anemometer. The device uses a nickel-iron wire of a known 
temperature coefficient of resistance and is heated at a constant current 
above the boiling point of water. In dry air, the sensing wire is balanced 
in a bridge circuit at a steady-state temperature between the electrical 
heat supplied and the heat removed by dry airflow. Droplets impinging the 
sensing wire evaporate, cool the wire, and thereby cause changes in wire 
resistance. The resulting imbalance in the bridge circuit is proportional 
to the liquid water content. Since larger drops tend to be sliced by the 
wire and not totally evaporated, LWC is underestimated when significant 
concentrations of large drops exist in the clouds. 
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The external probe houses two calibrated resistance wires which are 
mounted in a 1.1 inch diameter tube. The sensing wire is exposed perpendicu- 
lar to the airstream and bisects the tube while the second wire is mounted 
with its axis parallel to the airstream direction and at the base of the 
tube where it is not subject to water drop impingement. The base mounted 
wire serves to compensate for variations in air speed, altitude, and air 
temperature and is connected as part of the bridge circuit. At times, water 
or ice collected at the compensating wire causes changes and variations in 
the zero offset of the device and results in unreliable estimates of LWC. 
Characteristics of the J-W device are presented in Table 2-4. 


Table 2-4 

SPECIFICATIONS OF THE JOHNSON-WILLIAMS 
LIQUID WATER CONTENT INDICATOR MODEL LWH 


Liquid water content range 

Droplet diameter 
response range 

Time response 

Accuracy 

Precision 

Power 

Envi ronmental 


Mechanical 

Weight 


0-3 g/m 3 or 0-1 g/m 3 
<~40 ixm 

1 second at 100 mph; 

0.25 gm -3 
0.1 gm" 3 

27 VDC; 15 amperes 

Temperature: totally compensated 

Altitude: 0-30,000 ft 

Speed: 100-400 mph 

Cylindrical probe: 1.1 in diameter 

2.1 in long 

Probe with electronics: approximately 12 lbs 


2 . 5 Rosemount Model 102 Total Tempe r ature Sensor 

The Rosemount temperature sensor makes use of a recessed platinum re- 
sistance element and a flow duct which is designed to achieve a consistently 
high ram recovery and to reduce the effect of droplet impingement by inertial 
separation. Flow within the duct is controlled by boundary layer bleed ports 
to improve time response and a right-angle baffle for separation of water 
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droplets from the sensing element. The housing is equipped with de-icing 
heaters without appreciable Increase in sensing error. The sensing element 
is connected to a bridge circuit where a high gain operational amplifier con- 
verts changes in resistance to a voltage analog of temperature. The sensor 
specifications are listed in Table 2-5. 

Table 2-5 

ROSEMOUNT MODEL 102 TOTAL TEMPERATURE SENSOR SPECIFICATIONS 


Type 

Platinum resistance thermometer: 

Recovery factor 

0.99 

Range 

-50°C to +50°C 

Output 

0 to 5 v DC 

Resolution 

Q.Q2°C 

Linearity 

0.02°C 

Time constant 

2 seconds 

Power (heater) 

260 watts nominal 

Qualified to MIL- 

■P-277230 (ASG) 


2.6 Dew Point Hyqrometers (E.G. & G. Model 137-C3 and General Eastern 
MpdeQ ITT 

Both devices when used in conjunction with a water separator is capa- 
ble of measuring dew point from an aircraft at speeds up to 300 knots and at 
altitudes up to 25,000 Xft. The Instruments utilize the thermoelectric or 
Peltier cooling effect to cool a mirror to the saturation temperature or dew 
point. A focused light on the mirror is reflected onto a photodetector which 
in turn drives the thermoelectric cooler when changes in the reflected light 
are sensed. The mirror is cooled when a decrease of condensate occurs. The 
system stabilizes on and controls about a particular dew layer thickness. 

The measurement of the mirror temperature under stabilized conditions is thus 
taken as the dew frost point temperature. Table 2-6 presents specifications 
of the two sensors. 


* 
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Table 2-6 


Range 

Accuracy 

Response 

Power 


Range 

Accuracy 

Response 

Power 


SPECIFICATIONS OF DEW POINT HYDROMETER 


E. G.& G. 137-C3 


Control Unit 


Sensor 


+50°C to -50°C 

±0.5°C above 0°C 
±1.0°C below 0°C 

2°C/sec 

115 VAC 50/1000 cps 
20 watts 


+70°C to -55°C 
±0.5°C 


2°C/sec 


General Eastern Model 1011 


Sensor 


-75°C to +50°C 

±0.2°C above 0°C 

+0.4°C between 0°C to -40°C 

*1.0°C between -40°C to -75°C 

2°C/sec 

115 VAC 50/500 Hz 
50 watts 
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3 . 


DATA SET DESCRIPTION 

A total of 209 tests were conducted during the January 1981 testing 
program at the NASA LeRC IRT. For each test point, cloud measurements were 
taken in the IRT test section from a J-W device, a Leigh IDS, and two of the 
PMS probes. Eight additional test points obtained on June 8, 1981, supplement 
the data set. The testing in June was conducted to evaluate possible effects 
of the steam tube wrappings on the PMS probe measurements. The Leigh IDS was 
not available for these tests. Measurements were taken under a variety of 
tunnel conditions and included variations in tunnel airspeed, temperature, 
cloud LWC and cloud median volume diameter (MVD). Airspeeds ranged from 50 
to 285 mph while temperatures varied from about -20°F (-29°C) to 68°F (20°C). 
Based on early tunnel calibrations with rotating cylinders, the clouds produced 
'In the IRT were characterized by LWC and MVD ranging from about 0.5 to 5 g/m 3 
and 10 to 30 jum, respectively. Due to the large matrix of testing conditions 
and to time constraints of the program, each test condition was not measured 
by all of the available PMS probes. 

3.1 Data Recovery Percentages 

Table 3-1 presents a summary of the various PMS probe configurations 
in the IRT and data recovery information for each of the indicated cloud mea- 
surement devices. The test points are identified in chronological order for 
convenience. The eight test points (210-217) are listed separately since the 
two PMS probes and J-W device used in June were not identical to the probes 
used in January. The PMS probes are identified by type and a number from 1 
to 4 for reference purpose. In this report the PMS probes will not be identi- 
fied by organization since it Is not pertinent to the analyses. 

As shown, a large portion (~66 percent) of the J-W measurements from 
January were judged unacceptable and resulted from a loose wire connection 
in its control box which rendered the device as being either inoperative or 
intermittent. Measurements from the Leigh IDS were judged unacceptable if 
the device was not operating, operating intermittently, or operating beyond 
its full scale range. It is noted that seven of the tests were conducted at 
temperatures warmer than freezing on 21 January. Hence, since the Leigh 
device does not provide LWC estimates at these warm temperatures, these test 
points are not counted in the totals. 
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PMS PROBE CONFIGURATION AND DATA RECOVERY SUMMARIES 
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*The FSSP and ASSP used in June 1981 are identified as FSSP 1* and ASSP 2', respectively. These probes were not the 
same probes used in January 1981. 



Except for FSS'P 2, sufficient numbers of test points were recovered 
from the remaining PMS probes. Glass bead calibrations performed prior to 
each test day and review of the measurements revealed that FSSP 2 was not cali- 
brated to the manufacturer's specifications. Calibrations performed with four 
diffrent batches of bead sizes revealed that instead of peak counts being 
found at four different and widely spaced channel numbers, counts from two of 
the intermediate sized batches peacked at an identical channel number. Hence 
although three of the sized batches < 10-15 /xm, 25-35 ^m, and 35-45 /xm) had 
peak counts in their appropriate channel numbers (3, 8 and 13), the counts 
from the intermediate sized batch (15-25 /urn) peaked in channel 8 instead of 
the nominal channel 5, The consequence of this problem resulted in bimodal 
mass distributions with a mode always found in channel 3 and with the second 
mode ranging between channels 6 to 8. Measurements from the remaining PMS 
probes under identical tunnel conditions displayed single mode distributions. 

An example of a typical bimodal distribution from FSSP 2 is presented 
in Figure 3-1. For comparison purposes, an acceptable measurement from ASSP 2 
for the same test point and a rejected measurement from FSSP 1 are also in- 
cluded. The example of the rejected data point is presented since the exponen- 
tial form of the distribution is character! stic of practically all of the 
rejected measurements from the remaining ASSP and FSSP probes. For all of 
the rejected data from ASSP 4 on Janaury 19 (test points 141-209), the pre- 
test calibrations revealed significant changes in their measurement capability 
from the earlier tests. Whereas the counts from the bead lots of given size 
ranges peaked in their nominal probe channels early in the test program, they 
later fell in corresponding smaller numbered channels; thus indicating a 
broadening of channel band width and decrease in probe sensitivity. Cali- 
brations performed on both of these probes with sized glass beads between 
35-45 ^m revealed peak counts falling from probe channels 12-13 down to 
channels 8-9 during the later tests. The apparent broadening of channel band 
width of from 3 to about 4.5 ,um resulted in distributions of exponential from 
which were characterized by number or mass density decrease with channel num- 
ber (or particle size) increase. Although not verified, experience strongly 
suggests that probe deterioration in both these instances resulted from water 
being lodged in the optical lenses of both probes. 
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TEST POINT 177 (FSSP 1) 
TEST PGIUT 5 (ASSP 2) 
TEST POINT 5 (FSSP 2) 






The remaining rejected PMS probe measurements resulted from probe 
deterioration due to problems associated with their optical system. Rapid 
heating from about -4 to 29°F for the last test point (217) after the probes 
were "cold socked" in the IRT for several hours caused Tense "fogging" in 
both FSSP 1' and ASSP 2*. As for the remaining rejected samples, ice build 
up within the sampling tubes was observed to have obstructed the laser beam 
and thereby affected the measurements. These problems were readily detected 
by the operator in real time from noting the characteristic exponential dis- 
tributions on the display console. 

Subsequent review of the reduced data lead to rejecting an additional 
4 test points and was based upon noting that both the estimated LWC and MVD 
for a given tunnel test condition (given air speed and tunnel pressure set- 
tings) differed significantly from the other measurements taken with the same 
probe. For these rejected samples, the LWC estimates were lower by a factor 
of 2 to 3 while the MVD was also lower by about 4 to 6 urn. 

3.2 Tunnel Settings 

Cloud conditions (LWC and MVD) within the test section were controlled 
by separate air and water pressure settings of the spray system and by the air 
velocity through the enclosed and refrigerated IRT. The following empirical 
equations which were based on early work with rotating cylinders show these 
relationships. 


LWC(g/m 3 ) = 6.545^38.7 (P w - P a ) 


MVD (Mm) 


1.09(43,9 -VT748V) V38.7(P W - P a ) 

/ Jj- + 10 

\1.13 


+ 0.0052V + 4 


where 


P a = air pressure setting, psig 
P w = water pressure setting, psig 
V = tunnel air speed, mph 
20 < P a < 75 
P w < 110 

5 < (P w - Pa) < 65 
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Table 3-2 presents a summary of the IRT cloud conditions which shows 
the number of test points as functions of given air and water pressure set- 
tings and tunnel air speeds. As shown, the greatest majority (212 test 
points) of the tests were conducted at 24 primary pressure (fixed air and 
water) settings where a minimum of at least 2 test points were obtained. 

Five of the remaining tests were conducted at five different pressure set- 
tings (ID 25). In four of the fixed pressure settings (ID 15, 18, 19, and 
20), tunnel air velocity was also varied systematically. Hence it is indi- 
cated that probe testing was subject to 54 different cloud conditions for 
which each is listed information on the control parameters, P a , P w , and V. 

Temperatures are also tabulated and serve merely to present ranges 
over which each of the cloud conditions were run. Further information on 
the number of test points as function of temperature range is summarized in 

I 

Table 3-3. Here it is shown that the majority of tests were conducted at 
temperatures between -20 to 19°F and that the largest number (89) of test 
points were obtained at temperatures between 10 to 19°F. Had all of the 
measurements been judged acceptable, it becomes quite apparent that it was 
not possible to obtain measurements from each of the PMS probes for all 
test conditions (fixed cloud and temperature). 

3.3 Data Selection 

Measurements taken ove~ periods ranging from about 20 to 300 seconds 
were processed to obtain representative estimates of each test condition 
which was maintained in the IRT test section for durations ranging from 
about 2 to 15 minutes. The sampling periods were chosen to represent the 
early steady state cloud conditions as revealed by time plots of key cloud 
parameters (LWC and MVD). Possible measurement deterioration with time in- 
crease due to ice buildup on the probes was the primary reason for selecting 
data early (as opposed to late) in the total spray period. Steady state con- 
ditions were usually measured at the test section some 40 to 60 seconds after 
the nozzles were turned on and properly adjusted. The nozzles were turned 
on after cooling (or heating) the IRT to a desired temperature and after in- 
creasing the tunnel air speeds from idle velocity to a designated velocity. 
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^Includes the 8 test points taken on June 8, 1981. 
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Table 3-3 


NUMBER OF TEST POINTS AS FUNCTION 
OF TEMPERATURE RANGE 


Temperature 

Temperature 

Number of 

Index 

Range 

Test Points 

1 

-20 to -11 

51 

2 

f 

O 

ct- 

0 

1 

i— » 

23 

3 

0 to 10 

44 

4 

11 to 20 

89 

5 

20 to 32 

3* 

6 

>32 

7 


Total 217 


♦Points taken from June 8 only. 
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Figure 3-2 presents a timi; plot of measurements taken at the test 
stand during the conduct of a typical test (211). After initiating the spray 
nozzles (0 sec), the spray cloud is shown to arrive at about 20 seconds and 
to reach near steady state at about 50 to 55 seconds. After approximately 
90 seconds and for times to approximately 320 seconds into the run, the PMS 
probe estimates are shown to decrease while the J-W trace remains at a nearly 
constant value. Although the largest decrease appears in the LWC trace of 
the ASSP, estimates of the FSSP on other occasions decreased at faster rates. 
The spray cloud was turned off after a five minute (at 300 sec) testing 
period and was sensed approximately 25 seconds later ( 325 sec) at the test 
stand as evidenced by the rapid decline of the indicated probe traces. 

3.4 Data Reduction Procedures 

Procedures employed in the reduction of the PMS probe measurements 
consisted of initially converting the raw counts in each of the 15 size 
channels to a number density N-j by the following equation: 

c i 

Ni = A.j x V x EFF x At 

where 

Ci = raw counts recorded in each channel number, i, in a time 
interval , At. 

Aj = sample area for each channel number. 

V = true air speed. 

EFF = efficiency factor 

For the ASSP and FSSP devices, Ai remains constant for all 15 size 
channels, while the constant EFF is determined from the ratio of the sum of 
the raw counts in 15 size channels and the separately recorded total strobe 
counts (ST) during the time interval At. 


EFF 


15 



ST 
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For probes without the velocity reject feature, EFF is assumed to be equal to 
one. 

The constant EFF also maintains a value of 1 for data from the OAP. 
However, in these optical array devices, the sample area varies with channel 
number as: 

Ai = AD x (K - i) x L 

where 

AD = channel width. 

K = N-l, 

N = number of sensitive elements in the photodiode array. 

L = depth of field. 

For the larger channels (8-15), L is larger than the probe aperture 
distance of 6.1 cm and hence remains fixed to it. Depth of field compensation 
for the smaller channels (1-7) is thus computed from the following equation: 

A^ = A-j x i = 1»7 

where Ai to the right of the equal sign is computed from the previous equation. 

Once determined, the number density values provide means for computa- 
tion of the following bulk parameters from integration of the data. 

1. Total number density, Nj 


15 

N T = I Ni 
i = 1 


2. Liquid water content, LWC 


LWC 



15 

Z 

i = 1 


N i 


D i 


3 
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where 

D-j = midpoint size of each channel, Dj = i x AD 
P w = density of water. 


Median volume diameter, MVD 


D = MVD .15 . 

MVD = i N, D.- 3 = I IN, D, 3 
D = D i 2 i = 1 1 

Table 3-4 presents a summary of the sampling characteristics of the 
PMS probes tested in the IRT. The band widths listed were determined from 
the glass bead calibrations performed prior to each test. 


Table 3-4 

PMS PROBE SAMPLING CHARACTERISTICS 


Probe 

AD 

(Mm) 

Ai 

(mrn^) 

N 

Veloci t\ 
Reject 

FSSP 1 

3.0 

0.243 

_ 

No 

FSSP 2 

3.0 

0.254 

- 

No 

ASSP 1 

3.5 

0.672 

- 

Yes 

ASSP 2 

3.0 

0.811 

- 

Yes 

ASSP 3 

3.5 

0.459 

- 

No 

ASSP 4 

3.0 

0.614 

- 

No 

FSSP l 1 

3.9 

0.547 

- 

Yes 

ASSP 2 1 

3.3 

0.653 

- 

Yes 

OAP 1 

2.0 

eq 

22 

- 

OAP 2 

2.0 

eq 

22 

- 
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All cloud data were corrected for local velocity differences from 
bockage effects from the probe supports and traverse hardware in tunnel test 
section without the probes in place Although no significant velocity differ- 
ences were measured at the sampling points of the PMS probes, pitot-static 
pressure measurements disclosed overspeeding at the Leigh IDS and J-W hot 
wire by factors of 1.18 and 1.13, respectively, and have been appropriately 
handled in the data reduction software. 

3.5 Data Volume Description 

Data listings in the form of a drop size distribution report for each 
test condition and PMS probe measurements are presented in a separate data 
volume as part of the contract requirements. An example is presented in 
Figure 3-3. Table 3-5 presents a summary of the information listed and 
identifying keys for the computer generated listings. It is noted that the 
IRT cloud expectations were computed from the equations given in Section 3.2 
for which the independent variables are also listed. Air and water presses 
settings, temperature, and air speed for each test point were obtained from 
the tunnel operator's log sheet. In addition to the written comments, nega- 
tive values are also listed for some of the data judged unacceptable. For 
measurements from the J-W and Leigh devices, a default value of -0.88 g/nr* 
is listed if they were inoperative while a negative sign is assigned to the 
estimate if they functioned intermittently. Missing relative humidity infor 
mation for the test period is indicated by a default value of -99 percent. 

To complement the distribution reports, the data volume also includes 
a computer generated test point summary for each of the data tapes listed in 
Table 3-1. The summaries precede the distribution reports. Figure 3-4 
presents an example of a summary for Tape 3. Identifying keys for informa- 
tion summarized appear in Table 3-6. The test point summaries are also 
presented in Appendix A. 
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Figure 3-3. EXAMPLE OF DROPLET SIZE DISTRIBUTION REPORT 


Table 3-5 
KEY FOR TEST POINT 

ORIGINAL PAGE S3 
OF POOR QUALITY 

LISTING 

! 

1 

. if 

1 

* t 

Information 

Key 

Units 


Test point number 

TEST POINT 

_ 


Date 

DATE 

MO/DA/YR 


Run number 

RUN ID 

- 


Start time 

START TIME 

LST 


Sample duration 

DURATION 

sec 

Ii 

IRT Settings 
Air speed 

TAS 

mph 

1 

b 

Temperature 

TEM 

F 

if 

Relative humidity 

RH 

0 

ij 

Air pressure 

AIR 

psig 

ii 

Water pressure 

H20 

psig 

fr 

MVD 

MVD 

microns 


LWC 

LWC 

g/m 3 


LWC from J-W device 

JW 

g/m 3 


Coefficient of variation* 

CV 

% 


LWC from Leigh device 

LEIGH 

g/m 3 


Coefficient of variation 

CV 

% 


Probe identifier 

PROBE 

name 


Size distribution/per channel 

Mid-point diameter, Di 

DI 

microns 


Normalized size, Di /MVD 

DI/MVD 

- 


Number density, Ni 

NI 

/cm 3 


Number density per band width, Ng 

ND 

/cm 4 . 


Mass density, M^ 

Mass density per band width, Mg 

MI 

g/oK 

jj 

MD 

g/cnr 


Percent mass, Mi /LWC 

i 

Cumulative percent, mass, 2 Mi/LWC 

i = 1 

MI/LWC 

% 

Yf 

ii 

SMI /LWC 

% 

H 

n 

li 

Probe LWC, 2 Mi 

PROBE LWC 

g/m 3 

i { 

;i 

Coefficient of variation 

CV 

% 

Ii 

!! 

Probe MVD 

PROBE MVD 

microns 

! 

ii 

Coefficient of variation 

CV 

% 

|i 

Probe total number density, Ny 

NT 

/cm 3 

i| 

Coefficient of variation 

CV 

% 

ij 

ii 

Data quality comments 

COMMENTS 

- 

ij 

it 

1 

*The coefficient of variation is defined as 

the ratio of the standard devia- 

. 1 

tion to the mean of the parameters computed from the one-second 

data 

• 


recordings; CV = cr/mean. 
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Figure 3-4. EXAMPLE OF TEST POINT SUMMARY 


Table 3-6 



KEY FOR IRT TEST POINT SUMMARIES 


Information 

Key “ 

| Units 

Data tape 

TAPE 

- 

Run number 

RUN ID 

- 

Day 

DA 

- 

Start time 

TIME 

Local 

Sample duration 

DT 

sec 

Tunnel air speed 

TAS 

mph 

Tunnel temperature 

TEM 

F 

Tunnel relative humidity 

RH 

% 

Pa 

PAIR 

psig 

Pw 

PH20 

psig 

Tunnel LWC expectation 

T-Q 

g/m 3 

Tunnel MVD expectation 

DV 

microns 

J-W LWC* 

J-Q 

g/m 3 

Leigh LWC 
PMS probe 

L-q 

g/m 3 

Probe name 

PROBE 


LWC 

Q 

g/m 3 

MVD 

DV 

microns 

/cm 3 

n t 

NT 

Total raw counts/sec 

t RCT 

xlO- 2 

Gamma distribution coefficient, a 

ALP 

- 

Gamma distribution coefficient, \ 

LAM 

xl0 -3 /micron 

Average probe activity 

ACT 

% 

Average probe efficiency 

EFF 

% 

Data quality comments 

REMARKS 

- 

*The coefficient of variation (CV) for 
beneath the probe estimates. 

the cloud measurements 

are listed 
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4. RESULTS AND DISCUSSION 

Test point measurements judged acceptable serve as the basic data set 
for the analyses. The following subsections present the results of the follow- 
ing analyses tasks. 

* IRT cloud repeatability. 

* IRT cloud uniformity. 

* Effect of steam tube wraps on spectrometer measurements. 

* Temperature and velocity effects on spectrometer 

measurements. 

* Instrument intercomparisons. 


4.1 IRT Cloud Repeatability 

For obvious reasons, demonstration of cloud repeatability is essential 
to the remaining analyses tasks. For the purpose of obtaining quantitative 
estimates of cloud repeatability, the analyses evaluates the sample to sample 
variations with respect to the mean by the following standard statistical 
equations . 


x = 



s 


2 


n 

£ 

j = 1 


(Xj - x) 2 /(n - 1) 


CV = _5_ x 100% 

X 

where 

Xj = a sample variable 
x = sample average or mean 
n = number of samples 
s = estimate of the standard deviation 
CV = coefficient of variation 


37 


Hence low values computed for the coefficient of variation indicate high 
repeatability and vice versa. 

Table 4-1 presents information on averages and on the temporal varia- 
tion (CV) of the various measurements for given tests from data collected (1) 
during the time when the spray cloud passed through the test section (spray) 
and (2) during the intial time period of the spray (sample). Both of these 
corresponding data sets had identical starting times. The statistics were 
based upon samples computed for each second of recorded data and for the indi 
cated durations (n = number of seconds). For a given instrument and a given 
test, the following observations are made. 

• Based on the statistics given for the accretion (Leigh and 
J-W) devices, temporal variations in cloud LWC during a 
given spray remains small; typically between 4-8 percent. 

• Although the statistics from the accretion devices indicate 
higher absolute variations with increase in average cloud 
LWC amounts, their percentages are relatively constant and 
hence considered not affected by different airspeed or 
temperature settings in the IRT. 

• The sampling durations chosen are sufficient to provide pre- 
cise estimates of the average spray conditions as evidenced 
by the near identical values listed for the spray and sample 
averages of the accretion devices. Of the five test condi- 
tions presented, four of the sample estimates are within 
about 1 percent while the fifth (test point 61) is within 

7 percent of their corresponding spray average values. Based 
upon sampling considerations and the relatively small CV 
shown, it is quite apparent that the 20 to 50 second periods 
are more than sufficient to provide precise estimates of the 
average spray conditions. 

• The PMS probe estimates of the average spray LWC show signi- 
ficantly wider percent variations than those from the accre- 
tion devices and thus indicate that they are less reliable 
rather than that the cloud conditions varied widely during 
the course of the spray. Aside from the steady conditions 
indicated from the two accretion devices, support is pro- 
vided from the sample values which generally show higher 
estimated averages having significantly smaller CV values. 
These differences between the spray and sample statistics 
reflect decreases in measurement precision during the course 
of the spray which resulted from ice buildup within the 
sampling tubes as indeed were confirmed by inspection after 
test conclusion. 
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Based on the near identical CV values listed for the PMS 
sample estimates of cloud I.WC and those from the accretion 
devices, it is Indicated that the PMS probe measurements 
taken early during given tests are valid. Hence it is 
inferred from the sample estimates that the temporal varia- 
tions in cloud MVD were also small and for these test 
points in the range of from 1 to 4 percent. 


Measurement deterioration due to ice buildup in the PMS 
probes results in a larger percentage decrease in estimating 
LWC than in estimating MVD as is inferred by the correspond- 
ing sample and spray averages and is as expected since LWC 
is proportional to the cube of particle size. 


Table 4-2 presents additional information on the repeatability of the 
IRT spray clouds. Summarized information on both cloud LWC and MVD from 
measurements from each of the cloud instruments listed include: 


Table 4-2 

IRT REPEATABILITY SUMMARY FOR A GIVEN SPRAY CLOUD 


Instrument 

Variation (%) in LWC 

Variat 

:ion (%) i 

in MVD 

Riven 

Spray 

Given 

Day 

Test 

Period 

Yearly 

Given 
Spray ' 

Given 

Day 

Test 

Period 

Leigh 

4 

10 

12 



- 

- 

J-W 

6 

21 

15 

12* 

- 

- 

- 

FSSP 1 

15 

36 

19 

- 

3 

4 

4 

ASSP 1 

6 

20 

24 

- 

2 

3 

3 

ASSP 2 

10 

15 

- 

- 

3 

3 

- 

ASSP 3 

6 

20 

33 

- 

2 

7 

7 

ASSP 4 

9 

37 

- 

- 

3 

7 

- 


*0ne tunnel setting. 


1. Variation during a given spray is determined from the 
average of the CV values of all accepted test points. 

2. Variation during a given day which represents an evalua- 
tion of test point averages taken under identical tunnel 
conditions (fixed air and water pressure and airspeed 
settings) at different times on a given day. 
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3. Variation over the test period which represents an evalu- 
ation of test point averages taken under identical tunnel 
conditions on different days. 

4. Yearly variations which represent an evaluation of the 
year to year variations in the IRT spray clouds of a 
given set of tunnel settings. 


Information presented for items 2-4 above represents a normalized form 
of the coefficient of variation (NCV) which is based upon (1) normalizing the 
sample to sample variations of the estimates (xj = LWC or MVD) taken at a given 
tunnel condition, s, to their average, K s , and 2) pooling the normalized 
squared deviations of all possible given tunnel conditions, k, as follows: 


NCV = 



100 % 


n s 

2 Xj 
i = 1 


n s > 1 


For evaluating item 2 above: 

Xj = sample estimate taken at a given tunnel condition on a 
given day. 

n s = number of samples taken at a given tunne' condition on 
a given day. 



For evaluating the test period on daily variation: 

xj = average of all sample estimates taken at a given tunnel 
condition on a given day. 

n s = number of different days when samples were taken at a 
given tunnel condition. 
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For evaluating the yearly variation: 

Xj = average of all sample estimates taken at a given tunnel 
condition during a season. 

n s = number of different test periods (e.g., January and June 
1981) when samples were taken at a given tunnel condition. 


Based on this procedure, the samples were normalized to account for abso- 
lute cloud differences at the various tunnel settings and pooled to lend confi- 
dence to the variation estimates since in most cases with the PMS probes, only 
1 to 3 tests were conducted at a given tunnel condition over a two day period. 
Hence *11 of the accepted samples from a given instrument were not necessarily 
used in evaluating the effects listed under each column of the table. 

Prominent features of the summary table include the following 
observations : 


Larger LWC variations exist between sprays taken at fixed 
tunnel settings than during given sprays. Measurements 
from the accretion devices show average temporal deviations 
of 4 to 6 percent and normalized test point to test point 
variations of 10 to 20 percent. 

Generally larger LWC variations are indicated from the PMS 
probe measurements and infers that their estimates of LWC 
are less reliable than those of the accretion devices. 
Larger variations are also found between sprays of fixed 
tunnel conditions than during given sprays. 

Although slightly larger MVD variations are also shown to 
exist between sprays, the variations of between 3-7 percent 
are small and not significantly different than the averages 
of 2 to 3 percent shown for given sprays. For a MVD of 
25 ym, these percentages translate to absolute variations 
of less than about 1 yin. 

Cloud variations among tests conducted at given tunnel 
settings were similar irregardless of whether the tests 
were conducted on a given day or on different days as 
indicated by similar values listed beneath the given day 
and test period headings. 

Although a 12 percent yearly variation is listed, it is 
noted that the estimate is based only on one tunnel setting 
(k = 1) and on measurements taken with the J-W device in 
January and June 1981 (n s = 2) and therefore does not pre- 
sent a representative estimate of the yearly variation. 
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Based on these observations it is concluded that the IRT clouds are 
quite homogeneous (t 5 percent in LWC and ± 2-3 percent in MVD) during 
given tests and are highly reproducible (± 10-15 percent in LWC and ± 3-7 
percent in MVD). These estimates are considered not to be inconsistent with 
estimates provided by NASA which are based upon errors due to pressure and 
velocity settings of the IRT operating equations given in a previous section. 

NASA curves on drop size (MVD) error as functions of maximum errors of 
air and water pressure settings (t 2 PSI ) and of tunnel velocity (t H mph) 
are presented in Figures 4-1 to 4-3, respectively. Percentage error curves 
in LWC are presented in Figure 4-4 for maximum air and water pressure setting 
errors of i 2 PSI. Distinguishing features of the NASA error curves include 
the following: 


The drop size error due to an error in air pressure setting 
increases with decrease in air pressure and in tunnel 
velocity, and with increase in the absolute difference 
between the water and air pressure and amounts to less 
than + 2 fim for a ± 2 PSI setting error. 

The drop size error due to an error in the water pressure 
setting increases with decrease in air pressure setting 
and in tunnel velocity but with decrease in the absolute 
difference between the water and air pressures and amounts 
to less than ± 2 ixm for a ± 2 PSI setting error. 

The drop size error due to an error in the velocity setting 
increases with decrease in air pressure setting and in tun- 
nel velocity, and with increase in the absolute difference 
between the water and air pressure settings and amounts to 
less than ± 0.3 jiim for a t 2 mph setting error. 

For a + 2 PSI error in either the air or water pressure 
setting, the percentage error in LWC increases exponen- 
tially with decrease in the absolute difference between 
water and air pressure settings. For this setting error, 
a LWC error of t 5 percent is indicated for a 20 PSI 
pressure difference (PD) and increases exponentially to 
+ 20 percent at a PD of 5 PSI. The effect of a ± 2 mph 
velocity setting error has a negligible influence on cloud 
LWC and MVD. 


4.2 IRT Cloud Uniformity 

Cloud uniformity is concerned in this analysis task with cloud varia- 
tions along a direction perpendicular to the air flow in the IRT test section 
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ERROR IN DROP SIZE, MICRONS 



100 m.p.b. 


250 m.p.h. 
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DROP SIZE ERROR IN IRT DUE TO VELOCITY SETTING ERROR 
ERROR IN VELOCITY of + 2 M.P.H. (MAX ERROR) 





PERCENT ERROR IN L.W.C 



Ph 2 o - Pair, PSI 

Figure 4-4. LWC ERROR IN IRT DUE TO PRESSURE SETTING ERROR 
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and serves to determine whether or not cloud (LWC and MVD) differences exis- 
ted between the sampling locations of the two PMS probes. The probes were 
separated by approximately 8 inches. To attain this objective, measurements 
were taken at fixed and known distances from the center of the test section 
by moving the test stand along the cross-wind direction. For the evaluation, 
measurements made at fixed tunnel conditions on 8 January 1981 and on 30 Octo- 
ber 1979 (Boeing-Vertol tests) serve as the data base. Unfortunately, both 
the Leigh and J-W devices were inoperative during this segment of testing in 
January and were not installed during the Boeing-Vertol tests. 

Table 4-3 presents comparisons of the cloud measurements for these tests 
which were conducted at two different but fixed tunnel conditions. From the 
95 percent confidence interval (Cl) given for the five centered measurements 
and from the approximate 10 percent CV determined for the two off-centered 
estimates, it can be concluded statistically that the given cloud was not 
significantly different at + 12" from tunnel center on 8 January. Since 
only one test point measurement for a given tunnel condition was obtained at 
each of the test stand positions in the Boeing-Vertol tests, no statistical 

Table 4-3 

TUNNEL UNIFORMITY COMPARISONS OF 
CLOUD LWC AND MVD AT GIVEN DISTANCES FROM 
TUNNEL CENTER FOR FIXED TUNNEL CONDITIONS 


Date 

Test 

Stand 

Position 

No. of 
test 
points 

LWC 95% Cl 

( g/rn 3 ) (g/m 3 ) 

MVD 

(pm) 

95% Cl 
(pm) 



ASSP 2 

ESTIMATES 



1/8/81 

12" North 

1 

0.09 

21 



Center 

5 

0.17 0.10-0.25 

23 

19-27 


12" South 

1 

0.11 

19 

- 



BOEING VERTOL 

- ASSP ESTIMATES 



10/30/79 

24" North 

1 

0.60 

17 



Center 

1 

0.88 

19 

— 


24" South 

1 

0.87 

19 

- 
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conclusions can be given. However, based on the statistics given for the 
PMS probe measurements on cloud repeatability, it is strongly inferred that 
the LWC and MVD estimates at the given test stand positions were also not 
significantly different. These inferences are strongly substantiated by a 
large number of tests conducted by NASA IRT engineers during February 1981. 

In these independent tests, vertical bars (2 inch cylinders, 6 ft long) 
were placed at nine locations across the 9 ft width of the tunnel. Ice 
accretion was measured on these bars at a number of heights for several cloud 
and air speed settings. Since tunnel velocity remains essentially constant 
in the test section except for distances within 6 inches of the walls, cloud 
uniformity can be determined from contour plots of ice accretion. Based on 
these test results, NASA engineers concluded that the clouds remain uniform 
(±10%) in an approximate 2 ft high by 3 ft wide center section of the 
tunnel . 


4.3 Effect of Steam Tube Wraps 

The series of tests (210-217) conducted on 8 June 1981 were specifi- 
cally designed to determine whether or not the steam tube wraps (STW) had 
any significant effects on the measurement precision of the PMS scattering 
probes. A total of eight tests were conducted at a given cloud condition and 
included five test points with the wraps on and three without them. The 
tests were conducted at two tunnel temperatures; warm (23 to 29°F) and cold 
(-4°F). Test point 217 which was made without the wraps at the warm condition 
was rejected from the analysis. Subsequent rapid heating of the tunnel after 
an approximately two hours operation at the cold temperatures resulted in 
severe underestimating of particle sizes by both PMS probes due to fogging 
of the optics. 

Table 4-4 presents a summary of the measurements. For a given probe, 
the LWC estimates from measurements made without the STW are shown to be 
within the range of those made with the STW and therefore the small differ- 
ences (± 5 percent) in the averages shown are considered insignificant. In 
comparing the estimates of MVD, no difference is indicated in the FSSP mea- 
surements while only a one micron or 5 percent difference is indicates in the 
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ASSP measurements. Since this 5 percent difference in MVD is similar in mag- 
nitude to the test point to test point variations given in Table 4-2, it is 
also considered insignificant. It is therefore concluded that the STW did 
not affect the measurement accuracy of these scattering probes as expected. 

Table 4-4 

SUMMARY OF PMS PROBE MEASUREMENTS 
WITH AND WITHOUT STEAM TUBE WRAPS 


Test 
Poi nt 
Number 

Temper- 

ature 

(°F) 

Steam 

Tube 

Wrap 

J-W 

LWC 

(g/m 3 ) 

FSSP1 ' 
LWC 

(g/m 3 ) 

‘ ASSP2* 
LWC 
(g/m 3 ) 

FSSP1 T " 

MVD 

(pm) 

ASSP2' 

MVD 

(Pm) 

210 

29 

ON 

0.71 

0.71 

2.61 

22 

21 

211 

23 

ON 

0.77 

0.82 

2.39 

22 

20 

212 

-4 

ON 

0.77 

0.78 

2.19 

22 

21 

213 

-4 

ON 

0.80 

0.77 

2.26 

22 

21 

214 

-4 

ON 

0.83 

0.76 

2.17 

22 

21 

Average 



0.78 

0.77 

2.32 

22 

21 

215 

-4 

OFF 

0.77 

0.73 

2.48 

22 

22 

216 

-4 

OFF 

0.75 

0.72 

2.41 

22 

22 

Average 



0.76 

0.73 

2.45 

22 

22 

% Diff = 

/ 0N-0^F 
\ ON 

^ xlOO 

4 

5 

-5 

0 

-5 

Time plots 

of probe 

estimates 

of LWC and 

MVD for tests 

conducted at 


-4°F with and without the STW are presented in Figures 4-5 and 4-6, respec- 
tively. Although degradation is shown in the ASSP measurements in both 
tests, the adverse effect of not having the STW is shown (Figure 4-6) at 
approximately 240 seconds into the run when ice buildup practically disabled 
the probe. For the FSSP probe, no apparent change (compared to the J-W trace) 
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Figure 4-6. TIKE PLOT OF INDICATED PARAMETERS TAKEN AT -4°F WITHOUT STEAM TUBE WRAPS 



is indicated with time in the test without the STW. However, some deteriora- 
tion is shown in the measurements with the STW., The far greater deterioration 
in the ASSP was caused by interference to the laser beam from icing over of 
the water deflector protrusion which was located near the tip of the laser and 
withif/ its sampling tube. The FSSP devices do not have these obstructions in 
the sampling tube. To relieve this problem in some of the ASSP devices, the 
deflector pieces have since been removed and replaced functionally by warm 
air flow through a purge line terminating at the base of the sampling tube. 

As expected, since mass is proportional to the cube of particle size, 
probe icing is shown to result in greater underestimations of LWC than of MVD. 
In the extreme case shown of the ASSP measurements in Figure 4-6, a decrease 
of more than one order of magnitude is indicated (~270 sec) in the LWC esti- 
mate for a factor of two underestimation in MVD. For the FSSP estimates shown 
in Figure 4-5, the 2 pm or 10 percent change in MVD corresponds to a reduction 
in LWC of about 30 percent. Hence, during normal use and especially during 
moderate to severe natural icing conditions, it is indicated that the probe 
estimates of particle sizes or MVD would tend to be more precise than of the 
higher moment related parameters such as LWC. 

4.4 Effects of Temperature and Velocity on Probe Performance 

Of concern to users is knowledge of the measurement accuracy of the 
various cloud instruments under various operating and environmental conditions. 
The IRT provided opportunities to investigate two possible problem areas; 
namely the effects of temperature and velocity on instrument performance. 
Studies by NASA LeRC personnel suggest that all ice accretion instruments 
including the Leigh IDS systematically underestimated LWC at temperatures 
near the freezing point of water. This decrease in accuracy has been attrib- 
uted to the Ludlam Limit which is associated with a condition where less than 
100 percent of the supercooled deposits freeze upon impact as a result of 
heating caused by the release of the latent heat of fusion. Primary factors 
which determine the Ludlam Limit on the device are the ambient temperature, 
collection rate and perhaps droplet sizes. Since higher numbers of observa- 
tions were taken during the January program, it seemed appropriate to further 
examine this effect. Unfortunately, the tests were conducted at temperatures 
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which were either much warmer than the operating limit (<32°F) of the device 
or much colder (<17°F) than the temperature range of 20 to 32°F where instru- 
ment accuracy was suspect. 

For LJers with turbojet aircraft, information On the accuracy of the 
spectrometer measurements at airspeeds close to and beyond the upper limit 
set for by the manufacturer of about 250 Kts (~285 mph) is also desired. Of 
particular interest is to what degree are the particles undersized at or 
beyond the upper velocity limit. Unfortunately, tunnel blockage by the test 
stand and associated hardware limited top airspeeds to approximately 285 mph. 
Nevertheless it was still felt worthy of investigation. 

In order to evaluate both of these effects on measurement precision, 
absolute standards must first be known. In both of these investigations, the 
cloud LWC and MVD set forth by the tunnel operating equations (Section 3.0) 
were chosen as references. Since the operating equations indicate that the 
cloud LWC and MVD are also dependent on tunnel velocity, and since ony a few 
runs were made at a fixed cloud condition but with varying temperature, the 
ratios of the estimates to those given by the operating equations were strati- 
fied into given velocity and temperature ranges for evaluation. Table 4-5 
presents the stratification intervals used in the evaluation with accompanying 
index values. 


Table 4-5 

DATA STRATIFICATIONS 


Index 

~1WJV) 

TAS (MPH) 

1 

£-16 

50 

2 

-15 to -6 

100 

3 

-6 to 0 

150 

4 

0 to 10 

200 

5 

10 to 20 

225 to 27S 

6 

20 to 30 

>275 

7 

>30 

- 
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The results of this evaluation procedure are presented in Figure 4-7 
for measurements taken with the Leigh and J-W devices. The means and standard 
deviations (error bars) are plotted against the corresonding index values. 

For the Leigh measurements, no systematic change with temperature increase 
is shown in the mean ratio values as expected. Except for the observations 
made at a tunnel velocity of 50 mph, the remaining average ratio values dis- 
play no apparent dependency on velocity. At the lowest tunnel velocity, the 
Leigh device consistently underestimated the tunnel values. The random fluc- 
tuations of the means also indicate that the J-W estimates were not affected 
by the tunnel temperatures. Although a general increase in the values of the 
means is shown with speed increase, the wide error bars indicate that the 
effects of tunnel velocity were insignificant. These observations are con- 
sistent with the few runs made at fixed temperature and pressure settings but 
at varying tunnel velocities. Plots of these test point estimates for both 
these devices appear in Figure 4-8. 

In the MVD evaluations (Figure 4-9) of the PMS probe measurements, the 
means also remain fairly constant indicating that temperature changes did not 
bias these estimates. However, the means are shown to increase systematically 
with increases in tuitrel velocity. This dependency is further exemplified in 
Figure 4-10 from the test series conducted at fixed temperature and pressure 
settings and show that while the tunnel values based on the operating equa- 
tions gradually decrease with speed increase, the probe estimates remain 
fairly constant for speeds exceeding about 100 to 150 mph. Although these 
comparisons may perhaps suggest possible bias in the probe estimates, they 
may likewise infer possible errors in the operating equations. Discussions 
on these comparisons are presented in the following section. 

4.5 Instrument Intercomparisons 

Of primary interest to the total analyses effort is perhaps the com- 
parability of the measurements from the various instruments. Aside from 
establishing the compatability of the measurements from given clouds, the 
measurements also lend means to reevaluate the calibrations performed with 
the rotating cylinders. Table 4-6 presents regression analysis results (by 
the method of test squares) for the various instruments listed and includes 
the following information: 
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Table 4-6 


INSTRUMENT COMPARISONS 
LINEAR REGRESSION ANALYSIS SUMMARY 
y = a + bx 


Instruments 

n 

i : 

a 

b 

r 

e 

Compar- 

isons 

J-W vs Tunnel* 

80 

0.75 

0.24 

0.28 

0.26 

LWC 

Leigh vs Tunnel 

154 

0.39 

0.60 

0.79 

0.21 

LWC 

Leigh vs J-W 

57 

0.45 

0.61 

0.51 

0.25 

LWC 

PMS vs Tunnel 

121 

0.17 

0.55 

0.58 

0.72 

LWC 

FSSP 1 vs Tunnel 

37 

0.21 

0.09 

0.65 

0.09 

LWC 

ASSP 1 vs Tunnel 

22 

0.02 

1.28 

0.87 

0.37 

LWC 

ASSP 2 vs Tunnel 

10 

-0.01 

0.57 

0.95 

0.27 

LWC 

ASSP 3 vs Tunnel 

29 

0.34 

0.81 

0.87 

0.51 

LWC 

ASSP 4 vs Tunnel 

23 

-0.51 

0.81 

0.89 

0.21 

LWC 

PMS vs Ptf? 

121 

0.0 

1.00 

0.65 

0.66 

LWC 

PMS vs Tunnel 

121 

13.1 

0.50 

0.57 

2,6 

MVD 

FSSP 3. vs Tunnel 

37 

15.6 

0.47 

0.69 

1.8 

MVD 

ASSP 1 vs Tunnel 

22 

12.8 

0.39 

0.62 

1.8 

MVD 

ASSP 2 vs Tunnel 

10 

10.2 

0.60 

0.96 

0.9 

MVD 

ASSP 3 vs Tunnel 

29 

13.4 

0.50 

0.50 

3.3 

MVD 

ASSP 4 vs Tunnel 

23 

8.7 

0.70 

0.78 

1.8 

MVD 

PMS vs PMS 

121 

0.0 

1.00 

0.79 

1.9 

MVD 


*The tunnel estimates are based on the IRT operating equations. 
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1. n = number of data pairs. 

2. a = intercept of the linear equation of the form 
y = a + bx. 

3. b = slope of linear equations. 

4. r = correlation coefficient. 

5. e = standard deviation of regression. 

It is noted that for a given comparison, the estimates from the instru- 
ment listed first represent the dependent variables, while the estimates of 
the instrument listed second represent the independent variables. Hence, for 
the listed instruments, Leigh versus J-W, y and x are represented, respectively 
by the measurements from the Leigh and J-W devices. For these comparisons, 
the tunnel values were determined from the IRT operating equations. 

Comparisons with the Leigh and J-W devices included all of the accepted 
test point estimates. The PMS scattering probe comparisons, however, included 
only those measurements from given tunnel conditions (fixed airspeed and pres- 
sure settings) where the measurements from at least two different probes were 
deemed acceptable. This compromised data selection was used since comparisons 
between probes should be made of clouds having common microphysics properties 
and since only a few (four) sets of acceptable measurements from all five of 
the PMS probes were obtained under the same tunnel or cloud conditions. Based 
on this data selection criteria, measurements from 37 different tunnel condi- 
tions were compared. In addition, for each tunnel condition, the data from 
each probe is represented by the average of the corresponding test point 
estimates. 

Scatter diagrams comparing the LWC measurements from the Leigh and J-W 
devices against the tunnel estimates (IRT LWC) and against each other are pre- 
sented in Figures 4-11 to 4-13 and show that the measurements are generally 
in agreement with the calibrations performed with the rotating cylinders. As 
compared to the tunnel values, the measurements from both instruments are 
shown to fall on either side of the one-to-one line for IRT LWC of less than 
about 2 g/m^. Larger scatter is also evident of the J-W measurements in this 
range as is also indicated from the larger values of e listed in Table 4-6. 

For IRT LWC exceeding 2 g/m^, the measurements from both instruments fall below 
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Figure 4-11. 
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Figure 4-13. SCA 





the one-to-one line. For the Leigh measurements plotted at the IRT LWC of 
2.58 g/m 3 , the two lowest estimates were taken at a tunnel velocity of 50 mph 
where the device consistently underestimated the tunnel values. At the IRT 
LWC setting of 2.41 g/m 3 , the underestimations from the Leigh device is con- 
sidered to have been caused by incomplete droplet freezing. Although not con- 
sistent with the total data set, the decrease with temperature increase (from 
4 to 16°F) at this relatively high LWC condition suggests that the Ludlam 
Limit was reached at temperatures much colder than observed in the earlier 
program. This inference is not inconsistent with the observations made at 
the higher IRT LWC setting of 2.58 g/m 3 and the identical tunnel velocity of 
100 mph since the tunnel temperatures were much colder (-12°F). Studies are 
currently being actively pursued by NASA LeRC personnel on the heat balance 
of the Leigh device and should provide more definite answers on the various 
causal relationships associated with the Ludlam Limit. For the J-W measure- 
ments taken at these two higher LWC settings, it is found that the clouds were 
characterized by relatively large particles (MVD ~25 pirn) and, therefore, the 
device would be expected to underestimate the true LWC. It is also noted that 
the instrument performed rather intermittently (CV ~20%) at the 2.41 g/m 3 tun- 
nel setting and hence the test point estimate (or mean) would also be low. 

As compared to the tunnel estimates, the reduced scatter of the J-W 
points from the line of equality in the comparison diagram with the Leigh 
device suggests that part of the scatter results from the uncertainty of the 
operating equations. This observation is reflected in the higher correlation 
coefficient listed for the comparisons between the two instruments. It is 
further noted that the J-W measurements were generally lower than those of 
the Leigh device. 

The comparisons of the LWC estimates between the scattering probes and 
the tunnel are less encouraging as indicated by the wide range of values listed 
for the intercept and slope of the regression equation and by the corresponding 
scatter diagrams presented in Figures 4-14 to 4-18. The estimates from the 
individual probes either systematically overestimated or underestimated the 
the tunnel values and strongly suggests that in all but the estimates from 
FSSP 1, the discrepancies resulted primarily from uncertainties in their 
effective sampling areas. It is noted that at the two highest tunnel settings 
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Figure 4-15. SCATTER DIAGRAM OF ASSP 1 IWC VERSUS IRT LWC 
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(LWC >4 g/m 3 ), significant amounts of water existed in particles having sizes 
much larger than the upper size limit of these probes and therefore the probe 
estimates should be considered invalid. The weak dependency exhibited between 
the estimates of FSSP 1 and the tunnel may have perhaps resulted from a mal- 
function in its velocity reject circuit which was discovered during the pro- 
gram. When viewed (Figure 4-19) collectively, it is found that the individual 
probe estimates vary by large amounts. Although part of the scatter may be 
due to uncertainties in the absolute calibrations of the spray clouds, the 
variations between the individual probe estimates at given tunnel settings 
still remains large as shown in Figure 4-20. In this plot, the probe values 
are plotted against the means of the individual probe estimates for given 
tunnel settings. 

Although the data can be shown to be compatible to the tunnel values by 
applying appropriate correction factors, these results which were based upon 
use of the manufacturer 's specifications, clearly indicate that if the size- 
distribution measurements are used to estimate their bulk parameters (e.g., 

LWC), the probes should be calibrated against known standards. Also, due 
primarily to possible deterioration in measurement precision due to probe 
icing, it is further indicated that instruments such as the Leigh IDS or J-W 
device which make direct measurements of the integrated spectrum be used to 
complement the measurements. The relatively higher measurement precision 
offered by these accretion devices is indicated in Figures 4-21 and 4-22 
which show that the average test point values for the given tunnel settings 
fall on the one-to-one line with a minimal amount of scatter. Of these two 
comparisons, the average measurements from the Leigh IDS shows greater pre- 
cision and reliability. 

The average RVD estimates for the given tunnel settings for each of the 
probes are plotted against the corresponding tunnel estimates in Figure 4-23. 
Wide deviations appear in the tunnel size range of from about 10 to 20 jum 
where the measured estimates are shown to be significantly larger. A closer 
correspondence is exhibited at the larger tunnel sizes. These observations 
are reflected in the best fit values given for the coefficients of the indivi- 
dual probe comparison equations in Table 4-6 which shows large intercept values 
(8-15) and relatively small slope estimates of between 0.4 to 0.7. When the 
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Figure 4-19. SCATTER DIAGRAM OF INDICATED PROBE LWC VERSUS IRT LWC 
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Figure 4-20. SCATTER DIAGRAM OF INDICATED PROBE LWC VERSUS 
AVERAGE PROBE LWC FOR GIVEN TUNNEL SETTINGS 
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Figure 4-21. SCATTER DIAGRAM OF AVERAGE J-W LWC VERSUS 
IRT LWC FOR GIVEN TUNNEL SETTINGS 
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Figure 4-22. 5CATTER DIAGRAM OF AVERAGE LEIGH LWC VERSUS 
IRT LWC FOR GIVEN TUNNEL SETTINGS 







probe estimates are collectively compared against the tunnel estimates (PMS 
vs Tunnel), it is shown that not only does the correlation decrease but that 
the scatter (e) also increases significantly. 

Figure 4-24 presents comparisons of the probe estimates against their 
averages for given tunnel settings (PMS vs PMS) and shows that the estimates 
were in most cases (112 of 121 or 93 percent) within ± 3 /xm of each other. 

This error is consistent with the accuracy claimed by the manufacturer of the 
glass bead calibrations of ± one channel band width, which for these probes 
were between 3 to 3.5 pm. It is noted that part of LWC variations between the 
probe estimates can be attributed to these differences although not to the 
large extent shown in the LWC comparison plot against the tunnel values. Thus, 
although the measurement accuracy of the probes may be questioned from the 
comparisons presented in Figure 4-23, the good agreement found between the 
individual probes suggests that the MVD calibrations performed by the rotating 
cylinders may perhaps have been in error for some of the tunnel settings. 

Part of the disagreement between the probe and rotating cylinder estimates 
results from the fact that the probe MVD estimates did not decrease with tun- 
nel velocity increase as indicated by the operating equations. Unfortunately, 
the documentation of the procedures used in determining the MVD from the 
rotating cylinders was not available for this study and, therefore, further 
evaluations of the disagreement cannot be pursued. 

Figure 4-25 presents examples of measurements taken at a given tunnel 
setting by all five of the scattering probes. The volume or mass distributions 
have been modified to account for measurement differences in LWC and in MVD by 
plotting the percentage mass in each of the channels against the scaled parti- 
cle diameter, D/MVD. Although the Individual probe measurements have been 
normalized for these differences, it is shown that significant measurement 
differences can still exist. Whereas the distributions from probes ASSP 1, 
ASSP 2, and ASSP 3 can be considered as being comparable to each other, those 
from FSSP 1 and ASSP 4 are not. Hence it is strongly suggested that calibra- 
tion procedures be upgraded to also include means to evaluate the shape or 
form of the size distribution measurements. 
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Figure 4-25. NORMALIZED MASS DENSITY VERSUS SCALED PARTICLE 
DIAMETER DISTRIBUTIONS TAKEN AT A GIVEN TUNNEL 
SETTING BY INDICATED PROBES 
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4.6 Drop Size Distribution Form 

Of general Interest to this study was the determination of a distribu- 
tion function which best described the droplet size distribution measurements. 
For this purpose, the applicability of the gamma, log-normal, and Langmuir 
distribution functions were investigated. The approach consisted of evaluating 
the goodness of fit of the PMS probe measurements with the assumption that the 
given distribution functions were valid. The following presents brief descrip- 
tions of these functions and procedures used in their evaluation. 


4.6.1 Gamma Distribution Function (GD) 

The gamma distribution function is described by three parameters and is 
given as: 

N d = No exp (-XD) for D>0 

= 0 for D £ 0 


where 

Nd AD = number density between D and D + AD 
D = droplet diameter 

N 0 = displacement parameter, No > 0 
oi = curvature parameter, a > >-1 

X = scale parameter, X > 0 

D* = unit constant 


A desirable feature of the GD is the ease at which moment related bulk 
parameters of the distributions can through simple integration be resolved. 
Some of these parameters include: 

a. Total number density, Nj 

N t = r°° N D dD = N D“exp (- XD) dD 

Jo Jo 

- n r (o; + i) 
x a+1 

where r represents the gamma function and N for mathematical convenience 
represents N 0 /D*. 
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b. 


Total volume density, M 
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M 



N d D 3 dD 



D a+3 exp (-XD) dD 


= n r(g+) 

c. Liquid water content, LWC 

LWC - I p w M 
0 

d. Reflectivity factor, Z 


Z -jf" N^dD 



Another feature of the GD is that the MVD (or D 0 ) of a given droplet 
size distribution is uniquely defined by the curvature and scale coefficients 
as follows: 


MVD = D 0 = (a + 3.67)/\ (2) 

Since the GD is described by three parameters any three moment related 
bulk parameters of a given DSD may be used to estimate their values (Takeud' 
1978). However, this method commonly referred to as the moment technique may 
also be inefficient and therefore iteration techniques which minimize the 
variance of the data points about their estimates can also be used. 

The GD can also be normalized to account for sample to sample varia- 
tions in the values of the descriptive coefficients by the following algebraic 
manipulations to Eq. 1 and use of the indentity given by Eq. 2. 
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N exp (-XD) 


(3) 




No 

D a 


h_ 

D"N 


■ exp (-(a+3.67) D/Do) 


in — )/ (a+3.67) = -D/D 0 


D a N 


or 


Y = ~D/D 0 


(4) 

(4) 


It is indicated by Eq. 4 that if the DSD measurements (Nq) adhere to 
the GD, then plots on a linear graph between the unitless expressions given to 
the left of the equal sign and the scaled particle diameter D/D 0 will result 
in data points falling along a straight line having intercept and slope values 
of 0 and -1, respectively. 

For each DSD measurement, the descriptive coefficients of the GD were 
determined by the two different techniques mentioned earlier. Hence,, two 
separate evaluations on the applicability of the GD were performed. To im- 
prove upon the efficiency of the moment technique (Gamma I), four moment 
related parameters were employed and consisted of the following computational 
procedures. 

a = -b -V^~T 4ac (5) 

2 a 

where 

a = 1 - Q 

b = 11 - 5Q 

c = 30 - 6Q 

Q = LZ 

MF 

15 

L « £ N-j D-j 

i = 1 
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F 

Z 


15 

2 

i = 1 


"i D ? 


15 

2 

i = 1 


D 


4 

1 


2 N, D? 
i = 1 


X = [(q!+ 6) (a+5) (a+4) (a+3) (a+2) l/Z]l/5 


N 


_ M X 


a +4 


r (a+4 J 


( 6 ) 

(7) 


TH<' r a/» Ar\fl i«a f Ka A ( Cnmma T M a ■{ rf a/J a f aa wf a «m t aa 4 f a ati^ i a ap An f nr 
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linearized form of the GD given by Eq. 3. 


= N exp (-XD) 

[J 

or 

y = N'-XD (3a) 

where 

y = 1 n Ng - a 1 n D 
N' = In N 0 


Equation 3a serves as the basis from which simple linear regression 
analyses (method of least squares) resolves unbiased estimates of the descrip- 
tive coefficients. The iterative technique (see Takeuchi and Chein, 1V9 for 
further details) basically seeks a value of a for which the standard deviation 
of regression is minimized. Once minimized, values for N 0 and X are deter- 
mined from the "best fit" estimates of the intercept (N 0 = exp (N 1 ) and 
slope, respectively. In this study, a computed from Eq. 5 served as the 
initial value in the iteration process. The MVD is defined by Eq. 2 in both 
techniques. 
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4.6.2 Log-normal Distribution, LN 

If the volume (or mass) distributions are assumed to be log-normal, 
then the probability density function p(D) is given as: 



where 

p(D) • In = normalized volume density between D and D+AD, 

Dg = geometric mean diameter or antilog of average In D 

a = standard deviation of In D 

In order to objectively compare the goodness of fit of the measurements 
against that of the GD it is mandatory that Eq. 8 be linearized as follows: 

— — /in D-ln D Q \ 2 

In (crV27p(D ) = - ~ 1 ) 

\ o / 

-2 <y 2 In (<xV?7p(D)) = ^ln 

+ [-2 a 2 In (aV?7p(D))] 1 / 2 = InL 

D g 

or 

+ Y* = In JL 
D g 

exp (-Y 1 ) = D/D g for D/D g i 1 

and 

exp (Y') = D/D g for D/D g > 1 

Y = D/Dg = D/D 0 (9) 
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It is so indicated by Eq. 9 that if the volume distributions are in 
fact log-normal, the plots on a linear graph between the unitless expression 
given to the left of the equal sign and the scaled particle diameter D/D g 
will result in data points falling along a straight line having intercept and 
slope values of 0 and 1, respectively. Also, since the volume frequency dis- 
tribution is symmetrical about its geometric mean, the MVD is also identified 
by D g . 

It is further indicated that except for the difference in the sign of 
the slope, the values for both the intercept and slope and the expression, 
D/Dq, given to the right of the equal sign in both linearized forms of the 
GO and LN equations are identical. Hence, the goodness of fit of given mea- 
surements v.an be evaluated and objectively compared between the two assumed 
distribution functions by evaluating the root mean square (RMS) error of the 
data points from the straight lines. 

The coefficients of the LN and normalized volume density (Vj) were 
determined from each of the probe measurements by the following equations: 


Ni D? 

_T_ - 

2 Ni D 


(10) 


15 

1 n D q = 2 Vi In D-j 

y i = 1 


(ID 


a 



(V,. in o,. - In D g ) 2 



( 12 ) 


4.6.3 Langmui r Di stributions 

The Langmuir distributions are assumed normalized volume density dis- 
tributions about the scaled particle diameters D/D 0 as given (Langmuir and 
Blodgett, 1946) below: 


85 



Table 4-7 


°f P OOR buAU« 


FOUR ASSUMED DSD 


LWC 

(%) 

B 

D/D 0 

c 

(D/D 0 )l-5 

D 

(D/D 0 ) 2 -° 

E 

(D/D 0 ) 2 - 5 

5 

0.56 

0.42 

0.31 

0.23 

10 

0.72 

0.62 

0.52 

0.44 

20 

0.84 

0.77 

0.71 

0.65 

30 

1.00 

1.00 

1.00 

1.00 

20 

1.17 

1.26 

1.37 

1.48 

10 

1.32 

1.51 

1.74 

2.00 

5 

1.49 

1.81 

2.22 

2.71 


Common to these distributions is the given normalized volume density 
versus scaled midpoint diameter relationship of the B distribution. The 
remaining distributions are defined by different power values to which D/D 0 
of the B distribution is raised. Hence the Langmuir distributions are defined 
as a single parameter distribution function. 

Based on these considerations, the evaluation of the goodness of fit 
consisted of first determining the power values along given points of the 
normalized cumulative mass distribution of a given DSD measurement and second, 
computing the normalized RMS error of these values. The points along the 
normalized cumulative mass distribution were chosen to be at the midpoints 
as well as the intermediate points of the normalized volume distributions 
given in Table 4-7 and are listed in Table 4-8 along with the corresponding 
scaled diameters (DB) of the B distribution. Given D 0 , the determination 
of the power values at each of base 12 points (K) consisted of the following 
procedures: 

a. Determine through linear interpolation, the diameter, D K , 
at which the given percentage was found. 

b. Establish a scaled diameter, DS^ = Dj</D 0 „ 
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CUMULATIVE VOLUME DISTRIBUTION OF 
LANGMUIR B DISTRIBUTION 


Cumulative 

LWC 

(%) 

DB 

(D/D 0 ) 

Cumulative 

LWC 

(%) 

DB 

(D/D 0 ) 

2.5 

0.56 

65 

1.085 

5 

0.64 

75 

1.17 

10 

0.72 

85 

1.245 

15 

0.78 

90 

1.32 

25 

0.84 

95 

1.405 

35 

0.92 

97.5 

1.49 


c. From the corresponding scaled diameter (DBr) of the B 
distribution given in Table 4-8 determine the power 
value (E|<) by the following equation: 

log 10 (DS k ) (13) 

E k “ logJg(DB|J 

Given a number of DSD measurements, the normalized RMS error as defined by the 
following equation: 



Since the Langmuir distributions are characterized by a single param- 
eter, the deviations of the power values from their mean value at the given 
frequencies essentially provide a measure of the goodness of fit. In Eq. 14, 
the deviations are also normalized to their paramter estimates (or mean power 
values) as done in the RMS error determinations of the linearized froms of 
the gamms and lognormal distribution functions. Hence, the applicability of 
the three distribution functions can be objectively evaluated by comparing 
their RMS error estimates from a given set of DSD measurements. 
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4.6.4 


Normalized DSD 

Figures 4-26 to 4-32 present plots of the DSD measurements normalized 
to the three distribution functions. For the mere purpose of having the 
normalized measurements be inversely related to the scaled particle diameters 
as in the linearized GD plots, negative values of the normalized measurements 
(Y) of the LN are plotted against the scaled diameters. The normalized 
deviations of the power values of the Langmuir distributions are also plotted 
with respect to the normalized line of the GD and at the corresponding scaled 
diameters found for the given cumulative volume frequencies. 

Based on the least amount of scatter shown of the points about its 
straight line, it is indicated that the gamma distribution function best des- 
cribes the measurements. These observations are reflected in the P!MS errors 
listed in Table 4-9 for each of the assumed distribution functions. For the 
evaluations based on the entire spectra, it is shown that the RMS errors 
associated with the GD are lower than those of the LN by a factor of about 2. 
Also, of the two techniques used in estimating values for the descriptive 
coefficients of the GD, the estimates based on the iterative method resulted 
in slightly smaller errors as would be expected. In comparing the resulting 
errors between the three distribution functions for the mid-channel measure- 
ments, the average error resulting from the GD (Gamma I) approximation is 
shown to be lower than those associated with the Langmuir and LN assumptions 
by factors of 2 and 3, respectively. Although a slightly larger average RMS 
error is indicated from use of the Gamma II estimations of the descriptive 
coefficients, the difference is considered to be negligible. 

4.6.5 Comparisons with DSD from Natural Clouds 

Based on the close adherence shown of the measurements to the GD, it 
can be stated that the GD adequately describes the probe measurements taken 
in the IRT. Hence given DSD may be compared through their coefficient values 
provided that they all confirm to the GD. Figure 4-33 presents plots of 
measurements taken in natural clouds which have been normalized to the three 
assumed distributions and shows that these measurements are also best des- 
cribed by the GD. Characteristics of the ten DSD are summarized in Table 4-10 
which also includes evaluations of the goodness of fit. Since it has been 
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Figure 4-27. NORMALIZED DSD FROM ASSP I 
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Figure 4-28. NORMALIZED DSD FROM ASSP 2 
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Figure 4-29. NORMALIZED DSD FROM ASSP 3 
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Figure 4-30. NORMALIZED DSD FROM ASSP 4 
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Figure 4-32. NORMALIZED DSD FROM ASSP V 
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Table 4-9 

IRT CLOUD DSD FORM EVALUATION 




RMS Error (100% LWC) for 

Given DSD Form 

Probe ID 

No. 

of DSD 

Langmuir Log Normal 

Gamma I 

Gamma II 

FSSP 1 

83 

0. 230 

0. 117 

0. 098 

ASSP 1 

54 

0. 158 

0. 099 

0. 079 

ASSP 2 

32 

0.217 

0. 091 

0. 076 

ASSP 3 

48 

0. 217 

0. 126 

0. 078 

ASSP 4 

41 

0. 210 

0. 112 

0. 079 

FSSP 1' 

7 

0. 248 

0. 072 

0. 063 

ASSP 2* 

7 

0. 078 

0. 051 

0. 038 



Average 0. 194 

0. 095 

0. 073 



RMS Error 

(95% LWC) For Given DSD Form 

Probe ID 

Langmuir 

Log Normal 

Gamma I 

Gamma II 

FSSP 1 

0. 101 

0. 175 

0. 052 

0. 072 

ASSP 1 

0. 087 

0. 105 

0. 035 

0. 060 

ASSP 2 

0. 103 

0. 186 

0 o 061 

0. 066 

ASSP 3 

0. 156 

0. 200 

0. 075 

0. 073 

ASSP 4 

0. 127 

0. 179 

0. 079 

0. 076 

FSSP 1' 

0. 074 

0. 107 

0. 020 

0. C 4c 

ASSP 2' 

0. 059 

0. 055 

0. 012 

0. 028 

Average 

0. 101 

0. 144 

0. 048 

0. 060 


C- ~ ^ 
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Figure 4-33. NORMALIZED DSD OF NATURAL CLOUDS 
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shown that the estimates of LWC and MVD from the individual probe measurements 
differ under given cloud conditions, it is indeed difficult to compare the 
tunnel measurements against those made through the natural clouds. Neverthe- 
less, since absolute differences in the absolute mass distributions can be 
appropriately normalised and since the displacement coefficient, N 0 , is essen- 
tially removed in the volume frequency versus scaled particle diameter distri- 
butions, it may be appropriate to compare the values found for the curvature 
(a) and scale (X) parameters between these measurements. 

Table 4-11 presents a computer generated summary of the curvature (a) 
and scale (\) values of the GD for measurements taken at the various tunnel 
settings. It is shown that these DSD were characterized by oi values ranging 
from about 1 to 5 and by X values between 0.2 to 0. 5//jm. Also, for a given 
tunnel setting, variations in a and X between probes are indicated to be 
within about ±1 and ±0. 07/pm, respectively. In contrast, most of the DSD 
from the natural clouds shown in Table 4-10 had curvature values (Gamma II) 
of between 5 to 10 and X values ( X = — Hy j ^ ) between 0.5 to 1.0/^m. These 
differences are further exemplified in Table 4-12 wnich presents the mean 

i 

and standard deviations of the data sets. Where the IRT clouds show average 
a. values of aboaut 2.3 to 3.5, the natural clouds were, characterized by an 
of 6.8. It is noted that the range of a values from the natural cloud DSD 
are also consistent with the mean value of 8 found from approximately 500 
samples taken with oil slides from natural clouds in Russia (Khrgian, 1963). 
Figures 4-34 and 4-35 present examples of the measurements used in the com- 
parisons and basically show that except for Test Point 42, the normalized 
mass of the IRT clouds is distributed over a broader scaled size range. This 
observation is consistent with the f>ict that the IRT clouds where generally 
characterized by smaller curvature and scale values. In any event, the some- 
what broader irt drop size spectrum will not have a measurable effect upon 
the ice accretion, relative to the accretion in natural clouds. 
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Table 4-12 

SUMMARY STATISTICS OF MEAN 
AND STANDARD DEVIATIONS (o) 

OF CURVATURE AND SHAPE PARAMETERS 


Probe 

a 


Mean 

a 

! Mean 

a 

FSSP1 

3.4 

1.0 

0.33 

0.07 

ASSP1 

2.5 

0.9 

0,30 

0.06 

ASSP2 

2.3 

0.3 

0.25 

0.03 

ASSP3 

3.2 

0.8 

0.30 

0.06 

ASSP4 

2.4 

0.5 

0.27 

0.06 

Natural Clouds 

6.8 

1.7 

0.67 

0.21 
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MASS 



Figure 4-35. EXAMPLES OF NORMALIZED MASS DENSITY VERSUS SCALED 
PARTICLE DIAMETER DISTRIBUTIONS OF NATURAL CLOUDS 




5. CONCLUSIONS 

Tests conducted in the NASA LeRC IRT at different air speeds, pressure 
settings, and temperatures lead to the following conclusions: 

A. Instrument Comparisons: 

4 

It There was no apparent velocity effect to the PMS 
probe measurements at air speeds ranging from 
50 to 285 mph. 

2. The measurements from the LWC instruments (J-W and 
Leigh) agreed (little scatter) with each othe' and 
with the IRT calibrations. 

3. The five PMS probes gave an excessively wide scatter 
in estimates of LWC (0.2 to 2.0 compared to the IRT 
setting of 1 g/m 3 ), 

4. For given cloud conditions, the estimates of cloud 
MVD or drop size between the five PMS probes were 
within ±3 jim of each other in 93 percent of the cases. 

5. Significant failure rates were observed of these 
instruments especially wnen one considers that the 
instruments were required to be in proper operating 
condition for testing. Failure rates of 66 percent 
and 24 percent were observed in the 0-W and Leigh 
devices, respectively. Failure rates ranging from 
100 percent to 2 percent were observed at the six 
scattering probes installed in the tunnel. Some 

of the failures can be attributed to hardware prob- 
lems (e.g. v J-W and all of the measurements from the 
FSSP 2), while others were related to measurement 
determination due to probe icing and perhaps to 
wetting of the optics. 

B. Tunnel Concerns: 


1. The cloud conditions produced in the IRT are highly 
repeatable as indicated by the values of the co- 
efficient of variation determined for the following 
evaluations. 


Source of Variation 


CV LWC (%) CV MVD (%) 


During given sprays 
Between given sprays on given days 
Between given sprays among days 
Between given sprays among seasons 


±5 +2-3 

±10 ± 3-7 

±12 ± 3-7 

± 12 
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The PMS probe estimates of MVQ exceeded those of 
the tunnel values which were based on rotating 
cylinder calibrations. The "best fit" regression 
equation determined between the five PMS probe 
estimates (PMVO) and the tunnel values (TMVD) is 
given as: 


PMVD = 13.1 + 0.5 TMVD 


As compared to the log normal and Langmuir distri- 
bution functions, the gamma distribuation function 
best describes the drop size distribution measure- 
ments from the IRT and natural clouds. The drop 
size distributions from the IRT were slightly 
broader than those of the natural clouds as was 
also indicated by smaller values for both the 
curvature and scale parameters of the gamma 
distribution function. 



6. RECOMMENDATIONS 

The discrepancies shown between the measurements of the individual 
probes of given cloud conditions in the IRT strongly set forth the need to 
establish a test f aci 1 1 ty (or facilities) where Instruments can be calibrated 
against known cloud standards. Aside from providing means tv better define 
the sample area and band width of Individual probes, the resulting calibra- 
tions would add credibility to the measurements. An icing tunnel such as the 
IRT used in this study should be considered as a possible test facility in 
addition to those operated by the National Research Council (NRC) of Canada 
and by the Arnold Engineering Development Center (AEDC), Arnold AFS in 
Tennessee. 

Due to possible measurement deterioration from probe icing it is 
recommended that direct LWC measurement devices (e.g., 0-W and Leigh) such 
as used in this study complement the probe measurements during icing certi- 
fication flights. 

Based on the large differences found between the probe and tunnel 
estimates of cloud MVD, it is suggested that procedures used in estimating 
the MVD from the rotating cylinder measurements be reevaluated. Specific 
emphasis should be placed upon evaluating errors resulting from incorrect 
assumptions on DSD form. Also, the effects of air speed utilized in the 
present operating equations of the IRT should be reviewed in view of the 
discrepancies found in this study. 
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